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ABSTRACT 

In this work we investigate the accretion of cosmological fluids onto an intermediate- 
mass black hole at the centre of a globular cluster, focusing on the influence of the 
parent stellar system on the accretion flow. We show that the accretion of cosmic 
background radiation and the so-called dark energy onto an intermediate-mass black 
hole is negligible. On the other hand, if cold dark matter has a nonvanishing pressure, 
the accretion of dark matter is large enough to increase the black hole mass well beyond 
the present observed upper limits. We conclude that either intermediate-mass black 
holes do not exist, or dark matter does not exist, or it is not strictly coUisionless. In 
the latter case, we set a lower limit for the parameter of the cold dark matter equation 
of state. 
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1 INTRODUCTION 

Intermediate-mass black holes (IMBHs) are hypothetical 
compact objects with masses in the range 10^-10^ Mq, 
proposed to explain ultraluminous X-ray sources (ULXs), 
which exceed by 1-2 orders of magnitude the Eddington 
limit for the accretion luminosity of ste llar-mass black holes 
l|Angelini et al.l l200ll : iFabbiand [200^ . These objects are 
thought to form in the centers of globular clusters, where the 
stellar density is high enough to trig ger the runaway merg- 
ing o f stellar-mass compa c t remnants ([Portegies Zwart et al.l 
l2004l : iGiirkan et alll2004l : iFreitag et al.ll200a) . Careful mea- 
surements of stellar kinematics and density profiles of glob- 
ular clusters suggest indeed that some of these clusters may 
harb our IMBHs, although the evidence is still no t conclusive 
(e.g. lNovola et al.ll2008l : iBaumgardt et al.ll2005l ). The prop- 
erties of an extended central X-ray source in NGC 6388 also 
suggest an IMBH as a po ssible explanation for its central 
engine (|Nucita et al.l |2008| ) . The fact that many ULXs in 
nearby galaxies are associated with globular clusters pro- 
vides another argument for their existence. However, al- 
ternative explanations can be found for these observations, 
hence the existence of IMBHs is still a matter of discussion. 

The detectability of IMBHs depends on their masses, 
since the effects on their surroundings increase with mass. 
The mass of an IMBH is determined by its initial mass and 
the accretion of matter and energy from their surroundings. 
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This accretion can occur in discrete events, such as the in- 
fall of the stripped envelope of a star passing near the IMBH 
(e.g. Miocchi 2007), or by the continuous infall of any me- 
dia in which the IMBH is immersed, such as cosmological 
fluids (dark energy, dark matter, the cosmic microwave back- 
ground) or the intracluster medium of the globular cluster. 



The solution of the steady spherical accretion flow of 
a classical fl uid onto an isolated c o mpact object has been 
obtained by iBondi fc Hwi3 l|l944 l. iMichell (|l973) worked 
out the corresponding solution of the general relativistic 
fluid equations in the Schwarzchild metric , considering a 
fluid composed of massive particles, while iBabichev et al.l 
(2004j) developed the same model for massless particle 
fluids. In all cases, an important prediction of the models 
is that the accretion rate scales as the s quare of th e mass 
of the compact object. Applying Babichev et abl (|2004l ) 
model to an IMBH, one can readily see that the total 
mass change due to the accretion of dark energy during 
the entire life of the IMBH is negligible. However, this 
might not be the case for dark matter. Several authors 
have investigated the acc retion of dark matter b y su- 
permassive black h oles ([MacMillan fc Henriksen 20021: 



IZelnikov fc VasilevI 
■ 20051: 



2005allbl: I Munvaneza fc BiermannI 



Peirani fc de Freitas Pachecd 



20081 ). 



iGuzman fc Lora-Claviid ( 201lf ) computed time-dependent 
accretion models, finding that if dark matter is pressureless, 
it must contribute with a major fraction of the present 
black hole mass. The same reasoning might be applied to 
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IMBHs, suggesting that they might have accreted hundreds 
to thousands of solar masses in their hfetime. 

However, IMBHs are not isolated objects; they re- 
side in the center of massive stellar systems. The gravi- 
tational field of these systems might influence the accre- 
tion flow, enhancing the accreti on rate, an effect not taken 
into account in the mo d els of ISabichev et al.l l|2004 ) and 
iGuzman fc Lora-Clavii j (|201ll ). Hence, an investigation of 
the effects of the host stellar systems of IMBHs on the accre- 
tion flow seems relevant, as it would allow to make more pre- 
cise predictions about the accretion rates onto IMBHs and 
their final masses, which can be compared with the available 
observations. 

In this paper we explore simple models for the steady 
spherical accretion flow onto an IMBH located at the center 
of a globular cluster, taking into account the influence of 
the cluster. In Section [2] we develop our models, based on 
relativistic fluid dynamics in the gravitational field of both 
the black hole and the cluster. In Section [3] we apply them 
to compute the accretion rates for different kinds of fluids 
and we estimate the final masses of IMBHs. Finally, in Sec- 
tion |4] we discuss the impact of our results on the current 
knowledge of IMBHs and summarize the results of our work. 



2 THE MODEL 

Our model is based on the following hypotheses. The ac- 
cretor is an IMBH of mass M located at the center of a 
globular cluster. The mass distribution of the stellar sys- 
tem (cluster plus black hole) is assumed to be spherical, 
hence the space-time interval can be written as ds^ — 



■ ) , where r,9,4> are the usual 



spherical coordinates, t is the time and v and A are functions 
of r only. These functions are determined by the mass distri- 
bution of the stellar system, which in turn can be described 
by a single function m(r) that gives the mass m enclosed by a 
sphere of radius r. The Schwarzschild metric is the particular 
case in which m(r) = M, and hence v = —X = ln(I — 2A//r). 
Throughout this work, we use the natural system of units in 
which the gravitational constant G and the speed of light c 
are equal to unity. The accreted medium can be described as 
a relativistic fluid with an equation of state p = ujp, where p 
is its pressure and p its energy density. Its energy-momentum 
tensor is T^ji/ = {p + p)u^Uv — pg^iv, where g^i_v is the met- 
ric tensor and u'^ — dx'^ /ds is the fluid four- velocity with 
u'^u^ = f . 

We aim at computing the accretion rate onto the IMBH, 
M = — 47r limr->2A/ '■^T'o , for which we integrate the equa- 
tions of motion of the fluid, subject to the boundary con- 
ditions that far from the globular cluster the density poo 

and pressure Poo are constant, and the fluid is at rest 

I ' I t 

{uao = 0). Following iBabichev et al. (2004J, the integration 

of the time component of the energy-momentum conserva- 
tion law T.^l" = gives the energy conservation equation. 



{p + p){e " + "u^^ 



1/2 



(1) 



where Ci is a constant of integration. Another integral of 
motion is obtained by projecting the energy-momentum con- 
servation law on the four velocity, u^T.^" — 0. Integrating 
this equation gives the energy flux equation 



JS. 7^ e^^'^+^Kr'^ = C2, 



(2) 

where C2 is another constant. Combining eqns. [T] and [2] we 
obtain 



rp dp' 



p+p' I + e- 



V 2\ 

u) 



1/2 



Poo + Poo- 



(3) 



Considering the expression shown above, the accretion rate 
can be evaluated as 



M = -47r(poo -fPcx))C2 



(4) 



where we have used that in the limit r — )■ 2M the influence 
of the globular cluster vanishes and the metric approaches 
to that of Schwarzschild, for which X + v = 0. 

The constant C2 determines the accretion flux onto the 
black hole and can be ca lculated by fixing the parameters at 
any point. iMichell (f f 97fl ) has shown that the relativistic hy- 
drodynamical equations for accretion onto an isolated black 
hole have a critical point, which can be used to de termine 
C2- To obtain the critical point in our case we follow iMichell 
(I1971I ) and write eqns. [T] and [2] in their differential form, 
combining them into 



, A-i/ 2n— 1 A-iy 

+ e u ) e u 



-] du + 

ul 



-(e + e u ) (— e u + (\ 



A' ■iu' 2 ,^ , fv' A' 2 



ju e ) + 



dr = 0, 



(5) 



where we used the equation of state. We can see from this 
expression that if one of the bracketted factors vanishes, 
there is a turn-around point and the solutions are double- 
valued. Only solutions that pass through the critical point 
where both brackets vanish simultaneously are single valued, 
leading to 



i^'(r.)(f-a;)-^=0, 



(6) 



and 



Here, the subscript "c" stands for the corresponding vari- 
ables evaluated at the critical values. It is apparent from 
eqn.[7]that there is no critical point for uj < (dark energy). 
We will treat t his case s eparately later, and focus now on 
the case oj > 0. iMisner et al.. (,1970l) have shown that for a 
spherical stellar system. 



2m(r) 



r [r — 2m(r)] 
Hence is defined by the implicit equation 

m{rc) 2u) 



(8) 



(9) 



Tc I -I- 3aj 

Evaluating eqns. [2] and |3] at the critical point, we replace C2 
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Table 1. Model parameters. 



Cluster scale radius (ro) 

Cluster tidal radius (rt) 

Cluster concentration (cgc) 

Cluster galactocentric distance (Rgc) 

Black hole mass (M) 

Dark matter density at Rgc {p^^) 

Dark energy density (p^^) 



0.35 pc 
44 pc 

1.8 
3.1 kpc 
3000 Mq 
4.0 X 10-2ikg m" 
7.7 X 10-27kg m" 



in eqn. |4] to obtain the final expression in terms of known 
quantities, 



m (re) { 1 e2 



1~UJ 



(10) 



Note tiiat M is proportional to m^(rc) which, depending 
on the location of the critical point, can range from to 
(Mgc + M)^, where Mgc is the mass of the globular cluster. 
Hence, as Mgc 2> M, the accretion can be greatly enhanced 
if the critical point is located in the outer regions of the 
cluster. Eqn. [5] shows that for a given cluster, the location 
of the critical point depends solely on uj. 

In the next sections we apply this model to several fluids 
with diflferent values for the equation of state parameter cj, 
and analize the resulting accretion rates. 



3 ACCRETION RATES 
3.1 Perfect relativistic fluid 

As an example we consider the perfect relativistic fluid, for 
which UJ is of the order of unity and equal to the square of 
the sound speed in the medium in units of c. This model 
describes, for example, the cosmic microwave background 
(p — p/3). The description of the globular cluster plus black 
hole ma s s dist ribution m(r) was taken from the model of 
iMiocchil l|2007l ). with a scale radius ro = 0.35 pc, a con- 
centration Cgc = 1.8, and a central black hole mass o f 
M = 3OOOM0. These parameters, taken from lHarri3 (|l996l ). 
were chosen to represent the mass distribution of NGC 6388, 
one of the globular clusters with the strongest evidence in 
favour of the existence of an IMBH. 

Eqn. [9] shows the individual contribution of the black 
hole and the globular cluster in a very practical way. A plot 
of the left-hand-side of this equation for our model is shown 
in Fig.[T] The globular cluster contribution is the small peak 
at r/ro ~ 10, while that of the IMBH is the asymptotic 
growth at r — >■ 0. We can see that for relativistic fluids 
(lj 1) the right-hand-side of eqn.[9]takes values of the order 
of unity, several orders of magnitude higher than the maxi- 
mum globular cluster contribution. The critical radius then 
occurs at r ^ ro, where the black hole dominates the scene, 
and we can approximate the metric by that of Schwarzschild 
and solve for Vc, Uc and M . The problem in this case reduces 
to that of an isolated IMBH, a nd the accretion rate d epends 
only on , as predicted by iBabichev et al] l|2004l ). As a 
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Figure 1. Ratio m(r)/r for a typical globular cluster hosting 
an IMBH. The dashed and dotted lines define respectively the 
maximum and minimum values for 2a;/(l -|- 3lj) for which eqn. |9] 
has three solutions. 



toy example, for our model, the accretion rate of the cos- 
mic microwave background is M = 1.7 x 10~^^MQyr~^, 
a completely negligible value resulting from both the M- 
dependence of M and the extremely low photon density 
of the cosmic background during the matter-dominated era 
in which stellar systems for med. This value is in exce llent 
agreement with the results of lde Freitas Pachecol (120111 ). 



3.2 Cold dark matter 

Cold dark matter is a non-relativistic fluid with negligible 
thermal energy and pressure, hence p is essentially the rest 
energy density. For this fluid, the equation of state p = up 
with UJ can be used if we assume an isothermal flow. The 
same description applies to any perfect baryonic fluid at low 
temperature. Fig.[T]shows that for a fluid with a; < 4 x 10~*, 
the accretion flow is indeed modified by the presence of the 
globular cluster. 

For very low values of cj (lj < 10~^, below the dotted 
line in Fig. [T]), the critical point lies far from the cluster cen- 
ter (Fig. [2]), and m(rc) ~ M + A/gc. Hence, we expect the 
accretion rate to be greatly enhanced. There is also a tran- 
sition region in which eqn. [5] has three solutions (between 
the dotted and dashed lines in Fig.[l]). However, only two of 
them are physically admissible. If the middle one were the 
critical radius, the other two must be turn around points, vi- 
olating the boundary conditions on the flow velocity. Hence, 
only the inner and outer solutions are possible, but only one 
of them corresponds to the critical radius. This radius can be 
found by integrating the flow equation [5] for both solutions. 
Only in one case the integration converges, indicating that 
this is the correct value of rc. In the case of the cluster con- 
sidered here, the full integration of the flow shows that the 
critical point is the inner solution for uj > ujt = 1-5 x 10~^, 
and the outer one otherwise. At this transition value, the 
position of the critical point changes abruptly (Fig. [2]). This 
change, probably an artifact due to the simplicity of our 
model, does not affect our conclusions as we will show later. 

In Fig. [3] we show the results for the accretion rates de- 
rived from eqn. 1101 and using pac ~ p^^(? (see Table [T]), 
where is the dark matter density at the galactocen- 

tric distance of the globular cluste r, taken from the work of 
iKIvpin. Zhao fc Somervill^ l|2002l ). We can see the abrupt 
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• Inner solution 
■ Outer solution 
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Figure 2. Critical radius for a typical globular cluster/IMBH 
system, as a function of the equation of state parameter u). Circles 
and rectangles correspond to the inner and outer solutions of 
eqn. [9] respectively. 



change in the accretion rate, which foUows the correspond- 
ing change in rc. At cuu ~ 3 x lO"^'^, rc reaches the globular 
cluster tidal radius. For r S> rt the globular cluster cannot 
be considered an isolated object any more, as the tidal field 
of the Galaxy becomes dominant. The predictions of our 
model are not reliable in such a case. 

It is important to discuss the effect of the cluster (and 
thus the IMBH) motion relative to the d ark matter fiuid on 
the ac cretion rate previously calculated. iHovle fc LvttletonI 
l| 19391 ) have derived a formula for the accretion rate of a clas- 
sical fluid by an object moving through it, which is similar 
to that of Bondi &: Hovld (|l944l ) . According to the formula 
of iHovle fc LvttletonI (ligm " the relation between the ac- 
cretion rate M{v) onto an object moving with velocity v 
with respect to the fluid, and the accretion rate Mo onto an 
object at rest is 



M{v) = Mo 



3/2 ' 



(11) 



where Cg = is the sound speed in the medium. Com- 
puting the Newtonian limit for eqn. 1101 it can be easily seen 
th at in the limit a; — > the accretion rate approaches that 
of iBondi fc Hovld | 19441 ). Hence, we assume that the correc- 
tion factor given by eqn.[TT]can be used to derive M(v) from 
Mo in our models. In Fig. Owe show the accretion rate val- 
ues corresponding to cluster velocities v equal to 0, 100, 200 
and 500 km/s, which are typical for the motion of Galactic 
halo objects such as globular clusters. 

If we neglect the cosmological evolution of the dark mat- 
ter density, we can compute a rough estimate of the amount 
of matter AM accreted by the black hole along the glob- 
ular cluster lifetime M ~ 10 Gyr. For u < 10"^, M de- 
pends on m{rc, with rc ^ ro, hence the change in the ac- 
cretion rate due to the IMBH growth is negligible. Hence, 
we can compute AM ~ MAt > 10* Mq for any reasonable 
value of V. Thi s is no t consistent with pre s ent upper limits of 
iLanzoni et al] (|2007l '): iLutzgendorf et al] (|2007f ): ICseh et all 
l|2010t ) for the mass of the IMBH in NGC 6388. On the other 
hand, for uj > 10~® the accretion rate is low enough to avoid 
this discrepacy. This result is remarkable, as it implies that 
either the globular cluster does not host an IMBH, or there 
is a stringent lower limit on the cold dark matter equation 



Figure 3. Accretion rate for a typical globular cluster/IMBH 
system, as a function of the equation of state parameter ui, for 
different velocities of the cluster relative to the fluid. Dotted- 
vertical lines correspond to the upper and lower limits of the 
transition region. 



of state parameter, which cannot be null but at least of the 
order of 10~^. 



3.3 Dark energy 

In previous sections, we based our calculations on the ex- 
istence of a critical radius and, hence, a critical velocity 
for the flow. However, we can see from eqn. that for un- 
stable fluids (to < 0; se e ICarrofl. Hoffman fc Trodden! [2OO3I : 
iFabris fc MartiniliggTl for reviews on this subject) there is 
no critical point outside the blackhol e horiz o n (rc < 2M). 
Following Michel (1971) and Babichc v et all l|2004h . we as- 
sume that in this case the instabilities in the flow cause 
the growth of the fiuid velocity up to the speed of light at 
the black hole horizon, and use this as a boundary condi- 
tion instead of the critical values. With this assumption, 
th e accretion rate of dar k energy becomes identical to that 
of lBabichev et all (j2004l l. 



(12) 



This has a simple interpretation: in this regime, the pressure 
forces of a relativistic fiuid are of the order of the gavita- 
tional forces produced by the black hole near its horizon, 
hence much greater than those of the gravitational field of 
the cluster elsewhere, and the syste m behaves as if the blac k 
hole were isolated, as in the model of lBabichev et al.l \2004 ). 
Another interesting property of eqn. [T^] is that M is nega- 
tive when LJ < —1 (i.e., for phantom energy), leading to a 
decrease of the black hole mass. It is important to point out 
that this would be t he case only if the Generalized S econd 
Law is violated, as Ide Freitas Pacheco fc HorvarthI (|2007f ) 
have shown. For our cluster plus black hole model, and using 
the present dark energy density p^^c^ as p^x, (see Table [T|, 
we obtain 



M = 9.5 X 10" 



'M, 



M 



(13) 



As the coefficient of the right-hand-side of eqn. [13] shows, 
the accretion rate and the mass gained or lost by an IMBH 
of 10^-10* Mq during the past 10 Gyr is negligible. Hence, 
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we can conclude that the accretion of dark energy by an 
IMBH does not affect the evolution of the accretor mass. 



4 DISCUSSION AND CONCLUSIONS 

In this work we investigated the spherical, steady-state ac- 
cretion of cosmological fluids onto an intermediate-mass 
black hole at the centre of a globular cluster, taking into 
account the influence of the parent stellar system on the 
accretion flow. We also include in our models a correction 
for the effects of the motion of the black hole through the 
fluid. For relativistics perfect fluids (p p) and dark en- 
ergy, our results show that the presence of the cluster does 
not affect significantly the flow, hence our results coincide 
with those derived from mod els of accretion o nto is olated 
black holes, such as those of iBabichev et"al] (|2004l ). The 
conclusions of these authors hold, particularly that the ac- 
cretion rateJs_£ro£ortional to {poo + poo)M^ . As pointed 
out bv IBabichev et al] (|2004l ') this implies that, depending 
on the nature of the dark energy, the IMBH mass can ei- 
ther increase (for uj > —1) or decrease (for uj < —1). For 
uj = —1 the mass remains unchanged. For the particular 
case of IMBHs with masses in the range lO^-lO"* Mq and a 
typical globular cluster, our calculations show that the black 
hole mass is not affected by the accretion of either dark en- 
ergy or radiation from the cosmic microwave background. 

However, for the accretion of cold dark matter the sit- 
uation changes. We found that if dark matter is collision- 
less or has a very low speed of sound, the accretion rate no 
longer scales as the square of the black hole mass, but as 
the square of the mass inside the critical radius. This is due 
to the effect of the globular cluster gravitational field on the 
fiow. As the critical radius can be well outside the cluster 
core, this mass can be much greater than the IMBH mass, 
and hence the accretion rate can be enhanced by a factor 
of the order of 10*-10^. We also found that the accretion 
rate scales as c7^. These results lead to final IMBH masses 
greater than 10** Mq , well above the upper limits given by 
present observations. We point out that this enhancement is 
independent of the details of our simple model, particularly 
of the abrupt change in the critical radius observed as the 
speed of sound decreases. Even if this change were replaced 
by a smooth behaviour, the critical radius must approach 
the cluster outer regions as u decreases, greatly enhancing 
the accretion rate. Hence, there will still be a lower limit 
for UJ at which the mass accreted by the IMBH during its 
lifetime exceeds the observational upper bounds. Our mod- 
els estimate this limit at cj 10~^, which corresponds to a 
sound speed of Cs ~ 10 km/s. 

The last result has an important impact on our knowl- 
edge of IMBHs and dark matter. If IMBHs exist at the 
centres of globular clusters and dark matter has a low 
sound speed, IMBH masses should have grown to values 
beyond the observed upper limits. To restore the agree- 
ment with observations we are forced to assume that dark 
matter, if it exists, must have a sound speed at least of 
the or der of Cs ~ 10 km/s. Indeed. [Guzman fc Lora-Claviid 
arrive at the same conclusion investigating the ac- 
cretion onto supermassive black holes using more sophisti- 
cated, time dependent accreion models. Another possibility 



to check our conclusions is to investigate the accretion onto 
stellar-mass black holes. Some of these objects have known 
mas ses and their velociti es are well established. For exam- 
ple, iMirabel et al.l (|200ll ) calculated the orbit of the black 
hole X-ray nova XTE J1118-I-480. Computing the accretion 
rate for this black hole (with a mass of A/bh = 6.9Mq and 
a velocity of vbh ^ 140km/s, taken from these authors) we 
obtain negligible accretion rates. Hence, stellar-mass black 
holes do not seem to be useful for testing the dark matter 
equation of state. Thereby, if IMBHs were finally detected 
at the centres of globular clusters and detailed models of 
their accretion developed, these objects would turn into ex- 
traordinary tools to investigate this issue. 
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